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Abstract 

^ ■ The D meson production at forward rapidities in d+Au processes is calculated using 

0> . a pQCD based model, assuming that this treatment could be used as a baseline for 

distinct dynamical and medium effects. It is analysed how the nuclear effects in the 
QQ ' nuclear partonic distributions may affect this process at RHIC and LHC energies. 

Cp ■ An enhancement in the moderate qt region for RHIC, due to anti-shadowing in 

the nuclear medium, is found. Our prediction for LHC suggests that shadowing will 
suppress the D meson spectra for qx < 14 GeV. 

Q^' Key words: shadowing, nuclear collisions, small x QCD. 

^ ■ PACS: ll.80.La, 24.85.+p. 



o 



X 

^ I Relativistic heavy ion collisions provide an opportunity to study the QCD 

properties at energy densities several hundred times higher than the density 
of the atomic nuclei [1,2]. In these extreme conditions a deconfined state of 
quarks and gluons, the Quark Gluon Plasma (QGP), is expected to be formed 
in the early stage of the collision. These higher densities could induce a large 
amount of energy loss by gluon bremsstrahlung while hard partons propagate 
through the medium. The energy loss experienced by a fast parton may serve 
as a measure of the density of color charges of the medium it travels through 
[3,4]. In a dense medium, as the QGP, the energy loss may be huge. As large 
transverse momentum partons are produced very early in these processes, one 
expects that they can probe the early stage of the formed dense medium [5]. 
The mechanism of energy loss is thought to explain the observed suppression of 
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the high transverse momenta (qt) hadron spectra in central Au+Au coUisions 
at RHIC [6]. However, this feature of data could be explained [7] as well by 
saturation effects in the initial nuclear wavefunctions using the Color Glass 
Condensate (CGC) formalism [8]. 

In order to determine which mechanism is responsible for this suppression, 
d+Au collisions were studied at RHIC. At mid-rapidity the data [9,10,11] 
show an absence of jet quenching, which indicates that the observed high-^T 
suppression patterns in Au+Au collisions are not initial state effects encoded 
in the wavef unction of the beam nucleus, but are caused by final state in- 
teraction of hard partons with the produced dense medium. In order to test 
the consistency of the interpretation of these quenching effects due to energy 
loss in a deconfined medium, a comparative study of the attenuation patterns 
for massless and heavy partons was proposed [12,13,14]. The large mass of 
heavy quarks modifies the gluon bremsstrahlung, since it is suppressed for 
small angles 9 < ttiq/E [15], where tuq is the mass of heavy quark and E, its 
energy, implying different energy losses for heavy and light quarks propagat- 
ing in a dense medium. As a consequence, one can observe a softening in the 
light hadron spectra accompanying heavy quark jets, and a hardening of the 
leading charmed hadrons [12,13,16]. The production of heavy mesons are also 
affected by initial state effects, and their magnitude has to be estimated for a 
realistic prediction of energy loss of heavy quarks. 

In this work, the validity of the perturbative QCD and the coUinear factoriza- 
tion is assumed for RHIC kinematical regime, and this treatment is considered 
as a baseline to explicitate the presence of new dynamical effects in charmed 
meson production in relativistic heavy ion collisions. Focus is given to hadron 
(deuteron) - nucleus processes, since studies at this kind of interactions can 
provide important benchmarks for further measurements in nucleus-nucleus 
processes. In particular, the forward rapidity region is studied, where the nu- 
clear parton momentum fraction, ,X2, reaches smaller values and saturation ef- 
fects are expected to became important [17,18,19]. Since the x values reached, 
at RHIC energies, are not very small, it is necessary to keep in mind what 
are the predictions of the conventional QCD models which assume nuclear 
shadowing. In what follows our analysis concerns to the rapidity dependence 
of the nuclear modification ratios, it^AB, defined by 



^^^^^^^ - ABda^Jdyd^qr ' 

where y is the rapidity and A and B are atomic mass numbers, considering 
the EKS parameterization [20] of nuclear effects. Predictions for rapidity dis- 
tributions for D mesons at RHIC {^/s = 200 GeV) arc also calculated and 
the analysis is extended for LHC {^/s = 5.5 TeV). A comment related to 
the light hadron production is in order here. Currently, the description of 
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the experimental results in the central rapidity region can be obtained using 
a perturbative approach which includes the nuclear shadowing effects in the 
partonic distributions and an intrinsic transverse momentum of the colliding 
partons in order to reproduce the Cronin peak [21,22,23]. However, it is im- 
portant to emphasize that the existing conventional nuclear shadowing models 
cannot completely explain both Cronin effect and the suppression in the for- 
ward region, as well as the large value of that suppression. In contrast, in the 
framework of the Color Glass Condensate [8] both effects are predicted [17] 
to follow from the nonlinear evolution equation. If a similar scenario will be 
present in charmed hadron production is a subject of intense study [24,25,26], 
since the bulk of the heavy quark cross section comes from larger x values 
than for light hadron production in all rapidity range [27]. 

Hadrons with heavy quarks are a very important tool to study the properties 
of the strong interactions. Their large quark masses provide a scale which 
allows the use of perturbative QCD for computing production processes, since 
the long distance dynamics is effectively decoupled from the short distance 
dynamics. The value of the charm quark mass is in the limit of applicability 
of perturbative QCD, being a matter of discussion. In this framework, one 
can use the collinear factorization to calculate the heavy quark production. 
In leading twist, the semi-inclusive cross section factorizes into the product of 
gluon distributions, heavy quark fragmentation function and the hard partonic 
cross section. The application of factorization, however, is not evident for the 
kinematic regime where the heavy quark mass mq is much smaller than the 
center of mass energy, [26]. For instance, it has been claimed that in 
RHIC energy range, the x values reached by the nucleons are lower enough to 
justify the calculation with the nucleus assumed as a saturated dense partonic 
system (the CGC), with a characteristic saturation scale Qs, breaking the 
factorization of the process [24] . In particular, the heavy quark production, and 
consequently the D mesons production, was studied using this semi-classical 
approach in the Ref. [24,25], showing that the saturation phenomenon makes 
the spectra harder as compared to the PYTHIA prediction [24] . 

For jet production in a hadronic collision a pQCD-based model [22] is used 
in this work. In leading order pQCD, the pp' inclusive cross section (where p 
and p' stand for a proton (p) or a nucleon {N)) for production of a parton of 

flavour i = g,q,<l {q = u, d, s, . . .) with transverse momentum pt and rapidity 
y [28] is written as a sum of contributions of the cross sections coming from 
projectile (p) partons and from target {p') partons: 



daPP 



■iX 



da'P 



da'P 




yi,PT 



dyid^PT yi=y 
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where 




(3) 



(4) 



are interpreted, respectively, as the average flux of incoming partons of flavour 
i from the hadron p, and the cross section for the parton-hadron scattering. 
The rapidities of the i and j partons in the final state are labelled by yi and 1/2. 
In this model infrared regularization is performed by adding a small mass to 



the gluon propagator and defining mr — yPr'^Pi- The fractional momenta 
of the colliding partons i and j are xi,2 = ^(6=*=^* + 6=*=^^), with the integration 
region for 1/2 given by — log(yi/mT — e~^') < 1/2 < log(-\/s/mT — e^*). More 
details are given in Refs. [22,28]. 

For the charmed meson production at high energies the dominant subpro- 
cess is gg — > cc. This cross section da^^ /di can be found, e.g., in [29] and is 
proportional to q;s(//^), with ji — Qp — ^Jm\ + nig. The factor K in (2) is in- 
troduced in order to account for next-to-leading order (NLO) corrections and 
is, in general, energy and scale dependent [28]. For the parton distributions 
the CTEQ5 parameterization at leading order [30], evaluated at Qp, will be 
used and when the nuclear shadowing effects are considered in the calculation, 
the EKS parameterization [20] will be employed. 

Inclusive hadron production through independent fragmentation of the parton 
i into a hadron /i, is computed as a convolution of the partonic cross section 
(2) with a fragmentation function Di^h{z, Ql): 

^^pp'-feX daPP'^'^ . ^2^ 

where qr is the transverse momentum of the hadron h, yh its rapidity, and z 
the light-cone fractional momentum of the hadron and of its parent parton 
i. For details, see Eqs. (8)-(ll) of [28]. This pQCD-based model has been 
successful in describing the data for charged hadrons and neutral pions, at 
mid-rapidity [22]. In the low qx region, it was considered an intrinsic kx for 
the colliding partons, in order to correct the curvature of the hadron spectrum. 
However, since the interest here is the modifications due to nuclear shadowing, 
the intrinsic kx is not considered in this calculation. For the fragmentation 
function, the Peterson function [31] will be used with e = 0.043, as in Ref. [24] 
in the CGC framework. 
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Since in heavy quark production at high energies the dominant process is 
the gluon fusion, the cross section is strongly dependent of the behavior of 
the nuclear gluon distribution. Currently, there are several parameterizations 
in the literature which predict distinct behaviors and magnitude of the nu- 
clear effects in the gluon distribution and a recent comparison is given in Ref . 
[27]. For example, the EKS parameterization [20] has a strong anti-shadowing 
(i?^ = xG a/ AxGn > 1) at intermediate x {x ~ 0.1 — 0.2), due to momentum 
conservation constraint, and the EMC effect (i?^ < 1) at x ~ 0.2 — 0.8. For 
lower values of x, it presents shadowing (R^ < 1). On the other hand, the 
HKM one [32] , presents less shadowing at small x values and the EMC effect 
is not present at intermediate x. Furthermore, the momentum sum rule is 
underestimated by the fllJlNG parameterization [33], due to a strong gluon 
shadowing and a lack of anti-shadowing effect (For a recent NLO analysis 
see Ref. [34]). Due to these differences between the parameterizations, only 
bounds can be estimated for nuclear effects and the EKS one is used in order 
to provide a conservative estimate. The distinct effects in different x regions 
present in this parameterization create an asymmetry in the rapidity distribu- 
tion: at large negative rapidities, the nuclear momentum fraction, X2, is large, 
while Xi is small; conversely, the positive rapidities access small X2 and large 
Xi. Since our goal is to study the nuclear modifications, our analysis deals 
with positive rapidities. A similar study for light hadrons is presented in Ref. 
[35,23] (For an interesting discussion about this subject see Ref. [36]). 

In Fig. 1 we present the rapidity distributions of the D meson spectra at four 
distinct values of qt for d-l-Au collisions at y/s — 200 GeV. The dot-dashed 
curve, labelled CTEQ5, shows the prediction without nuclear shadowing and 
the solid curve, labelled EKS98, shows the predictions when the nuclear shad- 
owing is considered. An asymmetry is observed at low Qt, but disappears for 
higher qt- At qt = 2 GeV, a strong anti-shadowing enhances the spectra at 
negative rapidities, and shadowing suppresses it for positive ones. For increas- 
ing Qt, this asymmetry is weakened, with the X2 values at positive rapidities 
increasing, entering in the anti-shadowing region. At qx = 5 GeV, the rapid- 
ity symmetry is recovered. Higher values of Qt presents the reversion in the 
rapidity asymmetry, with enhancement of spectra at positive rapidities. 

At LHC energy {^/s = 5.5 TeV) and p-|-Pb coUisions, the cross section for the 
charm production probes the gluon distribution in the region of x > 3 x 10~^ 
for y < 3 [27]. In this region, the EKS parameterization, which is based on 
the DGLAP evolution equation and global fits of the DIS and Drell-Yan data 
above Q'^ — 1 GeV^ and x > 10~^, assumes that the nuclear gluon distribution 
behaves similarly to the nucleon one, which implies that the ratio keeps 
constant. Consequently, it does not consider any new dynamical effect associ- 
ated to the high density of the medium in this kinematic regime, which could 
modify xqa in comparison to xg^- This is a conservative assumption, since 
recent results for forward rapidity indicate that the inclusion of the saturation 
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effects is necessary. However, as our goal is to provide a baseline for future 
comparison, we use the EKS parameterization as input in our calculation. 
In Fig. 2 we present our predictions for the rapidity distributions for the D 
meson production at LHC energy. We have that at forward rapidities, the sup- 
pression in the spectra decreases with increasing g^, while the anti-shadowing 
dominates at large negative values. The crossover between the curves signal- 
izes the rapidity value where the nuclear shadowing begins to dominate and, 
with increasing g^^, this point gets closer the central rapidity. 

Our analysis regards to forward rapidities, since in this region the values 
reached for X2 becomes small enough to consider the saturation in the nu- 
clear wave function. The BRAHMS Collaboration has investigated the charged 
hadron production in d-|-Au collisions at forward rapidities, with the values 
7] = 1, 7] — 2.2 e rj — 3.2 [11]. Two of these values are considered to compute 
the evolution of the nuclear modification ratio Rab, defined in Eq. (1), in 
the transverse momentum. The prediction for mid-rapidity is also shown for 
comparison. Even thought the NLO corrections may affect the shape of the 
qt distributions, the higher-order corrections should largely cancel out in this 
ratio. Our results are presented in Fig. 3 for y/s = 200 GeV and for ^/s = 5.5 
TeV, where to compute the denominator d-l-Au processes were used in the 
first case and p-|-Pb processes in the later one. 

At RHIC, as shown in the left panel of Fig. 3, the spectra are enhanced in 
high-gr region, and the value of qt where the enhancement begins depends 
on the rapidity. If no effects are present, we expect the ratio as unity. At 
mid-rapidity, all the spectra is enhanced, weakening with increasing qt. With 
increasing rapidity, the spectra is suppressed at low Qt, but becomes enhanced 
for higher values of Qt, which is characteristic of the parameterization used. 
At fixed rapidity, x oc rriT/^/s ~ qx/^fs and the values of x increases with 
qr, entering the anti-shadowing region of the EKS parameterization. At very 
high r/i . wc expect that the ratio could fall below 1, due to EMC effect in the 
EKS. At RHIC, this result suggests that in the region where the validity of 
the perturbative treatment is expected, qx > 3 GeV, the D meson spectra will 
be enhanced due to nuclear anti-shadowing. This behavior is also present in 
charged pion production in same energies [23]. The preliminary open charm 
data from STAR Collaboration [37] at mid-rapidity show this feature in the 
region 1 CeV < < 4 GeV. In order to check our calculations, we have calcu- 
lated the D meson spectrum and verified that our results describe reasonably 
the experimental data [37] in the region of interest for this study (gr > 2 
GeV), underestimating the data in the region of lower qr, as expected, since 
we are not including an intrinsic transverse momentum. This result is not 
shown since our main focus is the nuclear modification factor RdAu, which 
can be described using a leading order calculation. On the other hand, in the 
calculation of the charmed hadron q^ spectra, the NLO corrections should be 
included, since it modifies the shape of the qr spectra as well as the normal- 
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ization of the cross section. Both effects largely cancel in the calculation of the 
ratio RdAu- It is important to emphasize that in the last Quark Matter Con- 
ference, the STAR collaboration has presented its preliminary results for the 
open charm spectrum in a broad transverse momentum region (0 < < H 
GeV) [38], with the measured open charm spectrum being much harder than 
the PYTHIA prediction, which is in line with the predictions from Ref. [24]. 
Such result indicate that the saturation phenomenon may be important for 
the heavy quark production at RHIC. However, more detailed studies related 
to the hadronization process are necessary before a definitive conclusion (See 
discussion in Ref. [38]). 

For LHC, as shown in the right panel of Fig. 3, the behavior of RpPb is similar 
for the three rapidities analyzed. At Iow-^t the spectra is suppressed, and the 
exact value where the enhancement takes place depends on the rapidity. At 
mid-rapidity, this point is qr ~ 14 GeV; for = 1, it happens at qt ~ 27 GeV, 
and this value increases for 77 = 2.2. Since the region to be studied at LHC is 
Qt < 14 GeV, this result suggests that a substantial suppression at positive 
rapidities is due to nuclear shadowing effects. Both panels in Fig 3 show a 
suppression in hadron spectra for qr < S GeV due to nuclear shadowing. 
However, the pQCD formulation might not be valid anymore for a quantitative 
calculation in this small Qt region. 

Finally, one expects that, because of their large mass, radiative energy loss for 
heavy quarks would be lower than for light quarks. It occurs due to combined 
mass effects [12,13]: the formation time of gluon radiation is reduced and their 
mass also suppresses gluon radiation amplitude at angles smaller than the ra- 
tio of the quark mass to its energy by destructive quantum interference [15] 
- the dead-cone effect. Due to these different energy losses, the ratio between 
hadrons with heavy quarks and with light quarks can provide a tool to investi- 
gate the medium formed in heavy ion collisions. The predicted consequence of 
this distinct energy losses is an enhancement of this ratio at moderately large 
transverse momentum, relative to that observed in the absence of energy loss 
(A recent analysis for LHC energies is given in Ref. [14]). 

In Fig. 4, results for the ratio between D mesons and tt, defined by 

nnAV~s,qT)^^, (6) 

are presented considering hadron-nucleus collisions, where the final state ef- 
fects, as energy loss, are minimal. The behavior of the ratio considering only 
the shadowing in the nuclear wavefunctions is presented. The ratio R^a^ de- 
fined as in Eq. (1), is calculated following Ref. [22] without the intrinsic trans- 
verse momentum. At mid-rapidity, for RHIC, the D meson production is more 
enhanced comparative to pions. This behavior also happens at large qt, for 
ri — 2.2. For LHC energies, the D meson production is smaller that the pions. 
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even at mid-rapidity. It suggests that the shadowing predicted in Fig. 3 at 
LHC for D mesons in p+Pb processes is stronger that the shadowing for pion 
production, due to the quadratic dependence on the gluon distribution present 
in the charm production. 

As a summary, the charmed meson production is studied using a perturba- 
tive approach. The rapidity distributions for d+Au processes at RHIC and for 
p+Pb processes at LHC are computed, and the disappearance of the asymme- 
try observed at low values of transverse momentum was found for increasing 
qt, since the X2 values increases with it. So, at high-g^ and RHIC energies, 
the D meson spectra is enhanced at forward rapidities. For LHC, the analy- 
sis predicts a suppression for positive rapidities due to nuclear shadowing, in 
the region qt < 14 GeV, even at mid-rapidity. We also studied the different 
behavior of charmed mesons and light hadrons in hA processes, where a min- 
imal energy loss are expected. Stronger nuclear effects for heavy quarks were 
found, which cause an enhancement for D mesons at mid-rapidity in RHIC, 
and their suppression for LHC, relative to pion production. This feature is 
based in a conservative perturbative approach, which assumes the validity of 
the collinear factorization and that the EKS parameterization is reasonable 
model for the nuclear effects. Although several points deserve more detailed 
studies, we believe that it can be used as a baseline for the CGC dynamics, 
expected to be present in this kinematical regime, as well as for future studies 
of jet quenching effects in AA collisions. 
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Fig. 1. Rapidity distributions for RHIC {^/s = 200 GeV) for distinct values of the 
transverse momentum. 
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Fig. 4. Ratio between D meson and pions, for RHIC and LHC energies and different 
rapidities. 
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